Switching the anomalous DC response of an AC-driven quantum many-body system 
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For a class of integrable quantum many-body systems, symmetric AC driving can generically 
produce a steady DC response. We show how such dynamical freezing can be switched off, not by 
forcing the system to follow the (arbitrarily fast) driving field, but rather through a much slower 
but complete oscillation of each individual mode of the system at a frequency of its own, with 
the slowest mode exhibiting a divergent period. This switching can be controlled in detail, its 
sharpness depending on a particular parameter of the Hamiltonian. The phenomenon has a robust 
manifestation even in the few-body limit, perhaps the most promising setting for realisation within 
existing frameworks. 



The coherent non-equilibrium dynamics of driven 
quantum systems offers a rich and largely unchartered 
field of many-body physics. While the physics of a pe- 
riodically driven single quantum particle has remained a 
lively topic of interesting research through past decades 
(see, e.g., QUI]), recent studies in periodically driven 
quantum many-body system have revealed surprising 
phenomena like dynamical many-body freezing [ME], 
appearance of slow solitary oscillations [7], superfluid- 
insulator and localization transition due to effective 
renormalization of the hopping [9HTT] which often dras- 
tically fail to conform to our classical intuition. Here we 
add a novel phenomenon to this collection, that of regime 
switching: at its most extreme, simply tuning a parame- 
ter in the Hamiltonian can lead to a transition between 
a completely frozen and an entirely dynamic regime. 

Our study contributes to the broader question of how 
quantum systems respond to non-adiabatic driving. A 
conceptually simple setting is provided by a quantum sys- 
tem driven periodically in time with a purely AC driving 
of frequency cj, with m z (t) denoting the "response" - 
a local quantity that varies in time under the influence 
of the driving. Let us assume [12] that the equilibrium 
value of this quantity vanishes when the driving field is 
zero, and that, like the driving, it is symmetric about 
zero over a driving period. Starting from the equilibrium 
state, purely AC driving must produce zero DC response 
over each cycle in the adiabatic limit (provided it exists, 
i.e. that there is no level crossing). However, if the driv- 
ing is non-adiabatic, the response can exhibit an average 
non-zero DC contribution even in the limit of infinite 
driving time. In an extreme case, quantum interference 
may dynamically stabilize the response around a nonzero 
DC value [6j [3 [14] . The (non-adiabatic) anomalous DC 
response is then defined as 



lim — 

t— )>oo T 



m z (t)dt . 



(i) 



Q can have strong non-monotonic behavior with a peak- 
valley profile in the space of driving parameters, fre- 
quency uo and amplitude ho [6[ [7] . For certain combina- 
tions of ho and u, the response remains frozen extremely 
close to its initial value for all time for any initial state. 



Here we discuss a very different but complementary 
phenomenon - how a non-zero Q can be made to disap- 
pear by simply tuning a coefficient in the Hamiltonian. 
This switching between Q finite and Q = is associ- 
ated with complete transfer of population between quasi- 
particle levels due to the driving, which surprisingly hap- 
pens on very long timescales even when the driving is 
arbitrarily fast. Under the special condition of maximal 
freezing, the transition approaches a discontinuous limit. 
The phenomenon has a dramatic precusor already in the 
limit of small system-sizes (> 6 spins). In fact, condi- 
tions for switching are less restrictive in this limit, and 
hence more easily observable in experiment. 

In the remainder of this paper, we first summarise 
the key observations for a quantum spin chain in a 
time-varying field. We then provide a detailed analytical 
theory, which agrees with a numerically exact treatment 
very well in the limit of fast driving. This is followed by 
an account of a finite-size version of the switching effect, 
together with a discussion of experimental realisability, 
and a concluding outlook. 



Regime switching: Consider the Hamiltonian 
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with <j a Pauli matrices, h z (t) = ho cos (cj£), and 
7 = Jx — Jy and S = J x + J y . 



(2) 



(3) 



We follow the transverse magnetization m z (t) = 
Z Si °f > taking the initial state to be fully polarized in 
the -fz-direction (m z = 1). Driving is strong and rapid, 



j, h > H,|J|. 
Fig. 0(a) 



shows m z (t) for a system of L = 10 4 
spins from a numerical solution of the time-dependent 
Schrodinger equation. After a transient, m z oscillates 
about a non-zero steady average, Q. This approaches 
zero as 5 — » 0. There is a concomitant substantial en- 
hancement of dynamics over long time-scales, visible in 
the increasing amplitude of oscillations about Q. 
The dependence of Q on driving parameters is shown in 
Fig. []Jb),(c). Besides the modulation of Q(oS) at fixed S 
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FIG. 1: Dynamical transition: (a): Exact numerical result for m z vs t for the quantum XY chain with uj — 10 over several 
periods for L — 10 4 . Inset shows the initial relaxation. We have taken ho = 20 and 7 = 0.5 (throughout the paper, unless 
otherwise noted), (b): 3d profile of Q vs S and uj. The sharp drop in Q at around S — (projected to the violet line on 
the base) and the peaks along the uj line (projected to bright yellow stripes on the base) are visible, (c): Q vs 5 for different 
uj for ho = 20. The figure shows sharp drop of Q at S = 0. The sharpness is a non-monotonic function of uj, approaching a 
discontinuous lim it for certain values of uj identified as the peak frequencies. The inset shows comparison between analytical 
formula (Eq. 10) for an infinite system with numerics for L = 100, averaged over > 10 4 cycles for uj = 10 (broader) and 
uj = U Q k + 10~ 2 = 7.256... (sharper). 



[6] [3 [14] , there is an even more pronounced dependence 
on S, with the response switched off, Q(uj, S = 0) = for 
vanishing S. The behaviour for small S is rich, exhibiting 
a tunable slope, the value of which diverges at the fre- 
quencies of maximal freezing uj p , yielding a discontinuity 
there: Q = at 5 = 0, and Q -)> 1 as 5 -> ± . 

Non-adiabatic evolution is of course not unexpected, 
since the gap closes at h z = ±7 in every cycle; moreover, 
uj, ho are much greater than the natural energy scales 
J x ,J y of the undriven system. However, disappearance 
of Q and enhancement of long-time dynamics as 5 — >• 
is surprising, particularly given the fact that reducing 
5 reduces the minimum gap, the smallest gap given by 
\ho cos cut + S\. 

Analytical theory: We reduce Hamiltonian in Q to 
the direct sum of 2 x 2 Hamiltonians in momentum space 
using Jordan- Wigner and Fourier transformations [15 : 

H = Y,i>\H k {t)i>k\ H k = [h(t)+df k }a z +A k a++A* k a- 

k 

(4) 

where k is the quasi- momentum, ipk = (c|(fc), C2(fc)) is 
the spinor representing the fermionic operators Ci^(k) in 
fc-space, a ± = (a x ±ia y )/2, S and h(t) (independent of 
k) are the parameters in the many-body Hamiltonian, 
and Afc = 7 sin k and Sfk = S cos k are functions of k. In 
the Schrodinger picture represents the Hamiltonian 
for the time-dependent wave function, which in turn is a 
direct product of /c-space wave functions: 

\m) = II \Mt)) = l[[Mt)\0k) + vk(t)\i k )]. (5) 

k k 

were {|0fc), |ljfe)} are the basis states in which the matrix 
form of Hk is expressed in Q. 

The Hamiltonian in Q represents a class of systems in- 
cluding well known integrable ID quantum models, such 



as the transverse Ising chain, ID p-wave superconduc- 
tor p~6j [17], or quantum spin ladders [18]. For instance, 
the Hamiltonian of a ID p-wave superconducting wire of 
spin-less fermions can be described by mapping 5 onto the 
hopping, 7 onto the pairing gap, and h(t) onto the chem- 
ical potential, with m z = 2nj — 1 related [16] the average 
fermionic occupation number, rif. A similar mapping is 
applicable to the Kitaev spin ladder in [T8] . 

The momentum space wave-function ([5| satisfies the 
following time-dependent Schrodinger equation 

i d ^ = H k \mh (6) 

and the time-dependent transverse magnetization reads: 

m z ( t ) = m)wt\m) = ^ E (M*)i 2 - V2) (7) 

k>0 

An exact analytical solution of Eq. [6] is not available, 
but for large uj, analytic solution is possible in the rotat- 
ing wave approximation (RWA) [6] [22] . Starting from a 
fully polarized (m z = 1) initial state, corresponding to 
Vk(t = 0) = 1 (a good approximation to the ground state 
of the Hamiltonian at t = for ho ^> 7, 5), 

K(*)| 2 = 1 - Mt)| 2 = J o(2WAf sin2 {<pkt) _ (8) 

where 

<t>k = ^Jl(2ho/uj)/\l + 5 2 cos 2 ^ . (9) 

For periodic boundary conditions, we have wave vectors 
of the form k = (2n + 1)tt/L, where n denotes positive 
integers. Averaging over time and integrating over k in 
the continuum (L — >• 00) limit, 

= 1*1 fio) 
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This formula (Fig.[T|)) is one of our central results. It very 
accurately reproduces the numerical results (Fig. J^:), in 
particular the switching off of Q for S = unless the 
Bessel function simultaneously vanishes. 
(Dis)appearance of Q: From Eq. (J8|, the max- 
imum amplitude of excitation for any given mode is 
j {2ho/uj)A k ^ w j 1 ^ c ] :l i s atmost unity, and in general 

smaller. This is what leads to a finite Q in our case, 
as is apparent from the expression of m z in Eq. (|7|). 
Note here the role of the apparent "DC" part 8 cos k in 
the 2x2 Hamiltonian (Eq. [6| , which plays the same role 
as would have been played by an additional DC part in 
the transverse field (albeit without the ^-dependence) . It 
is essential for the anomalous dc response: for 6 = 0, we 
have Q = for all modes and for all uj, however large! 

It is important to note that Q = emphatically does 
not reflect any adiabatic dynamics: eqs. (8j9) encode the 
drastic enhancement of dynamics of each /c-mode at 5 = 
0. At this point each mode has full oscillation with am- 
plitude \uk\ 2 (but not frequency fa) independent of both 
k and parameters of the Hamiltonian. Thus, the charac- 
teristic time of the full population transfer for a mode of 
quasi-momentum k, = 2ix jfa. Indeed, from Eq. ([9| it 
follows that T — >• oo for 5 = and cos k — >• 0, so that some 
modes oscillate on diverging timescales, leading to an 

asymptotic m z (t) ~ J (|2J (^) \t) ~ |^o (^) 
This form evidences the singular nature of the points 
(uj = uj p ,S = 0), where uj p are the zeroes of the Bessel 
function, as we discuss below. 
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FIG. 2: Short time vs long time dynamics: Discrepancy 
between small short-time and large long-time excursions for 
the dynamical case, 5 = 0, while, for 5 = 0.1, a fast initial 
decay is followed by long-time stability around the average 
Q (numerical result for L = 100). 



By contrast, the short-time behaviour offers little clues 
into what happens at long times (precluding a simple pic- 
ture based on naive repetition of the single sweep result) 
and in fact makes the dynamical curves look more static 
than the frozen ones, in keeping with the observation that 
switching involves times much longer than the sweeping 
period: Fig. [2] compares dynamics for S = and S = 0.1 



close to the sharp transition point uj 



10" 



For 



0, while the long-time dynamics takes the system 



far away from the initial state, at the end of first few 
cycles, the system returns remarkably close to the initial 
state. But for S = 0.1, the system remains quite close 
to the initial state for all later time, while its deviation 
from the initial state at the end of few initial cycles is 
systematically much larger than that for the S = case. 
Sharpness of the transition and the discontinuous 
limit: As seen from Figjl](b), (c), the sharpness of the 
transition varies non-monotonically with uj. For discrete 
values of uj, the limit is discontinuous: Q = 1 for S ^ 
and Q = for 5 = 0. As can be seen from Eq. (J8|, 
if 5 7^ 0, Q = 1 is maximal at the frequencies, where 
Jo(2ho/uj) = 0, and, up to remarkable accuracy all the 
/c-modes remain completely frozen for all time [6]. How- 
ever, if 5 = 0, all modes undergo full oscillation. Hence 
the limits Jo(2ho/uj)Ak —> and also 5cosk —> do not 
commute: we either obtain Q = or Q > 1, switching 
between complete dynamics and complete freezing. The 
discontinuity is of course derived under RWA, whose ac- 
curacy increases with uj, and it becomes exact as uj — > oo. 

The tunable sharpness of the switching might be in- 
teresting from the point of view of sensitive detectors of 
very small static magnetic field. As discussed in the pre- 
vious section, a small DC component of the transverse 
field plays the same role as a non-vanishing S. Now, if 
one sets (5 = 0, and tunes the system close to an ujf then 
a DC transverse field will result in a jump from Q = to 
Q = 1. Of course, the resolution time over which coher- 
ence needs to be maintained will grow as one approaches 
uj p '. Since the transition persists in the limit of small 
system consisting of few spins (see below and the Sup- 
plementary Material), maintaining coherence over the pe- 
riod necessary for detecting even very weak fields might 
not be unrealistic within the framework of cold atoms in 
optical lattices or ion traps (see, e.g. [20, 21]). 
Experimental possibilities: finite-size versions 
Various kinds of quantum spin Hamiltonians and their 
coherent dynamics are realizable in cold atomic systems 
in optical lattice (see. e.g., |20j El]). One important 
feature of the switching phenomenon that makes it par- 
ticularly convenient for experimental realization is its 
clear signature even at very small system sizes of just 
a few components (with periodic boundary condition). 
Suppose we are away from the transition condition (i.e. 
S ^ 0), and drive the system with uj, ho ^> 7, 5 close to 
the large system perfect freezing condition: some small 
systems, L = 3, 4, 5, 7, 8, 9 also freeze. However, at L = 6 
we see pronounced dynamics in the system, resulting in 
a strong dip in Q. Similar pronounced dynamics are ob- 
served at system sizes 10, 14 ... 2(2n + 1) with gradually 
decreasing amplitude (Fig. [3]). 

The origin of such strong dynamics even for 5 can 
be gleaned from Eqs. <Js]) and For L = 2(2n + 1) 
we have a pair of /c-modes for which Scosk = 0: again 
the term involving the Bessel factor responsible for the 
freezing cancels out for that particular pair of modes and 
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FIG. 3: Manifestation of dynamical transition at small 
system-size (L): The figure (main + inset (a)) shows (nu- 
merical results) drastic contrast of response dynamics for 
L — 5 and 7 with that for L = 6 for 5 = 7 = 1, uj — cj^ + 10 -2 . 
Even for 5 7^ 0, for small L strong dynamics is observed for 
certain system-sizes (L = 2(2n + 1), n — 1,2,3,...) repre- 
sented by the corresponding dips in Q in inset (b). These are 
the system sizes for which a mode k has cos/c = 0, which will 
have full oscillation independent of the value of uj (Eqs. [8j 
[9]) . The blue dotted line is Q = 1 — r/L with r is determined 
using the numerical value of Q for L = 6 (the first dip). 



account of how fast driving yields the appearance of long 
timescales non-monotonically. Given the applicability of 
the basic Hamiltonian studied here to a wide range of 
systems, and the appearance of switching even for small 
systems, experimental realisation does not appear an un- 
realistic prospect in the near future. 

Several further questions obviously spring to mind. 
Firstly, what are the necessary conditions for such be- 
haviour to occur beyond the soluble RWA - here, e.g., 
the point 5 = is special in that it exhibits an enhanced 
XY-like symmetry, as well as a degenerate Floquet spec- 
trum, and appears like a critical point with a diverging 
timescale - but it is not clear which of these fundamen- 
tally underpin such switching mechanisms in general. As 
a next step, progress towards studying the long-time dy- 
namics of non-integrable systems would be of interest to 
determine the stability of this phenomenology, this being 
hampered by the present inability to access the long-time 
dynamics of such systems. Finally, given the growing zoo 
of such nonequilibrium quantum phenomena, a system- 
atic zoology based on some yet to be identified organising 
principles is clearly a most desirable goal in the long term. 



they execute full oscillation. For small L, there are only 
few modes in total, and hence the contribution from the 
fully oscillating modes is substantial resulting in the dips 
in Q. As L increases, Q at the dips scales as Q(L) = 
1 — r/L, where r is a constant (inset of Fig. [3]). 

Since this small system-size signature of the transition 
does not require 6 — > 0, it can also be observed even in 
the transverse Ising chain, which is a special of the XY 
chain with (5 = 7 = 1. Indeed, shown in the figure are 
results for the Ising case where in fact no dynamics is 
actually expected in the large L limit, since in this case 
we have 7 = S = 1. We emphasize that the transverse 
Ising chain with time- varying transverse field has already 
been realized experimentally [23, 24 j. Different initial 
conditions: Freezing at uu p is maximal as there, every 
mode freezes individually [6], while for S = 0, every mode 
becomes dynamical. In that sense, freezing and switching 
are independent of the starting state, which only affects 
the amplitudes Uk(t — 0) and Vk(t = 0), which can be 
accounted for by a time-independent unitary transforma- 
tion. Such a transformation preserves Q = 0, although 
full population transfer would not occur in the modified 
basis. Thus switching and the concomitant discontinu- 
ous limit are essentially unaffected but the fully polarized 
initial state we have chosen provides a large Q away from 
the transition point, making the transition most clearly 
visible. 

Conclusion and outlook: We have presented switch- 
ing between (completely) frozen and dynamical non- 
equilibrium regimes under driving, with regime change 
effected by simply tuning a parameter in the Hamilto- 
nian. This set of phenomena is particularly intriguing on 
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